Orbital-quenching-induced magnetism in Ba2NaOsO( 
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PACS 71.27+a - Strongly correlated electron systems; heavy fermions 

PACS 71.70.Ej - Spin-orbit coupling, Zeeman and Stark splitting, Jahn- Teller effect 
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Abstract. - The double perovskite Ba2NaOs06 with heptavalent Os (d^) is observed to remain 
in the ideal cubic structure (i.e. without orbital ordering) despite single occupation of the t2g 
orbitals, even in the ferromagnetically ordered phase below 6.8 K. Analysis based on the ab initio 
dispersion expressed in terms of an Os f2g-based Wannier function picture, spin-orbit coupling, 
Hund's coupling, and strong Coulomb repulsion shows that the magnetic OsOe cluster is near 
a moment-less condition due to spin and orbital compensation. Quenching (hybridization) then 
drives the emergence of the small moment. This compensation, unprecedented in transition metals, 
arises in a unified picture that accounts for the observed Mott insulating behavior. 



Orbital physics in transition metal oxides has attracted 
a good deal of attention for three decades, with much of 
the activity focused on coupling to spin, charge, and lat- 
tice degrees of freedom in systems. The configura- 
tion is typically found in the early Sd transition metals 
(TMs), i.e. titanates and vanadates, and the associated 
phenomena - often revolving around the tendency to or- 
bitally order - are remarkably rich. The d^ configuration 
can also occur in the mid- to late-Sd TM ions, which are 
distinguished also by unusually high formal valence states. 

Unless fully itinerant, partially filled d shells lead to 
non-spherical ions that often are accommodated by orbital 
ordering (00), that is, ordered alignment of the filled or- 
bitals [1,2] in the manner that minimizes strain, electronic, 
and magnetic energies. [3] An outcome is that 00 has been 
identified as the driving mechanism in symmetry-breaking 
structural and magnetic transitions, [2] and may also be 
coupled to magnetism and charge order, the Mott insulat- 
ing d^ perovskite YTiOa being a well studied example. [4] 
Very different behavior is shown by the d^ system LaTiOs 
which has a structural distortion although evidence of 00 
has been difficult to obtain. A prominent explanation is 
that orbital fluctuations dominate, leading to a disordered 
orbital liquid ground state. [5,6] 

Double perovskite structure Ba2NaOs06 (BNOO) is a 
rare example of a heptavalent osmium compound, also rare 
because it is a ferromagnetic insulator [7-9] (Tc = 6.8 K), 
and it shows other perplexing behavior. Although it has 
a single 5d electron in the t2g complex that orders mag- 



netically, it shows no evidence of orbital order that would 
destroy its cubic symmetry. On the other hand, the sis- 
ter compound La2NaOs06 with high-spin Os configu- 
ration with nominally cubic symmetry is observed to be 
highly distorted, [10] which is ascribed to misfit arising 
from the small cation radius. The question of spin order- 
ing surely is a delicate one, since isostructural and isova- 
lent Ba2LiOs06 (BLOO) orders aniiferromagnetically in 
spite of a very similar Curie- Weiss susceptibility [7] and 
nearly identical volume. 

The question of formal valence, and identification of sev- 
eral material constants, can be identified from first princi- 
ples local density approximation (LDA) calculations using 
two all-electron full-potential electronic methods of FPLO 
and Wien2k. [11-14] The Fermi level Ep lies within the t2g 
bands, confirming the heptavalent nature of the Os ion. A 
crystal (ligand) field splitting of almost 5 eV separates the 
centroids of the Cg and t2g bands, reflecting the unusually 
strong 5d— 2p hybridization, and a gap of roughly 1.5 eV 
separates the t2g symmetry bands from the more tightly 
bound O 2p-dominated OsOg cluster orbitals, with gaps 
separating various degenerate molecular orbitals. [15, 16] 

Several complications must be taken into account to ob- 
tain a realistic picture of the electronic structure, (i) The 
strong hybridization of Os 5d orbitals with O 2p states, 
which has ramifications in this double perovskite struc- 
ture beyond the usual ones in oxides, (ii) Strong spin- 
orbit coupling (SOC) on the Os site has important con- 
sequences, (iii) Hund's (exchange) coupling effectively se- 
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Fig. 1: (color online) Nonmagnetic Os t2g band structure in 
LDA, in three Wannier orbital representation (tight binding 
[TB]), and LDA+SOC, relative to the Fermi level Ep. In LDA, 
the top of the t2g manifold is very flat near the P point and 
along the X-W line, resulting in very sharp peaks just above 
Ef in the density of states (see fig. 2). From LDA+SOC, the 
spin-orbit coupling constant ^ is 0.30 eV. 

lects the occupied orbital, (iv) Intra-cluster Coulomb re- 
pulsion, which is the obvious suspect for a magnetic insu- 
lating state in an open shell system. Finally, it is likely 
that quantum fluctuations may also play an important 
role in determining the observed behavior. While the in- 
terplay of spin and orbital physics has been addressed in 
considerable detail, [17] the large SOC in BNOO puts this 
compound in a distinct class relative to the early 3d sys- 
tems. 

The unusual aspects of BNOO then are: a Mott insulat- 
ing state based on a localized (f ion that nonetheless re- 
tains a cubic environment, and a small ordered moment in 
a 3D system where fluctuation effects typically are minor. 
It is most instructive to consider the various aspects of 
the electronic structure consecutively. We build an under- 
standing of the electronic structure and competing mecha- 
nisms in a different manner from Xiang and Whangbo [18] 
in their study of BNOO. 

Isolated t2g bands: cluster orbitals. The bands, with 
bandwidth W=l.Qb eV, form a self-contained system for 
the purpose of understanding the low energy, low temper- 
ature physics. The three bands can be represented accu- 
rately in a Wannier function basis with three first neigh- 
bor d-d hopping parameters (in meV) t^ = — 121,i7r = 
64, ts = 24 meV, with on-site energy e = 202. (For BLOO 
tc — —131 is the only change.) Note that ts is necessary 
to give the excellent representation of the bands shown 
in fig. 1. The corresponding total and atom-decomposed 
densities of states (DOS) are pictured in fig. 2, which re- 
veals the isolation of the narrow t2g bands. As pointed 
out by Erickson et al. [8] the overlap takes place between 
O 2p7r orbitals on neighboring OsOe octahedra; these O 
orbitals comprise the 'tails' of the t2g symmetry cluster 
orbitals. The Fermi level lies 0.5 eV (one electron) up 
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Fig. 2: (color online) Total and atom-projected densities of 
states (DOS) per formula unit for nonmagnetic state in LDA. 
The vertical dashed line indicates the Fermi energy. Inset: 
Blowup near Ef. It displays that contribution of six oxygens 
is nearly similar with Os character in the t2cj manifold. As a 
result, it is anticipated that the oxygens significantly contribute 
to magnetic moment. 

into the bands. The DOS shown in fig. 2 is highly non- 
symmetric about their center, with 2/3 of the weight lying 
in the upper half of the bands. 

As analyzed in some detail by Singh and coworkers 
[15, 16] in the context of pentavalent Ru in Sr2YRu06 
which also contains a closed shell ion (Y'^"'") in one of the 
cation sites, the OsOg unit should be approached as an 
isolated cluster. The Os t2g states generalize to cluster 
orbitals with the Os ^dxy state (say) bonded to the 
orbitals lying in the x — y plane. We calculate that half of 
the density of the t2g bands lies on Os, the other half dis- 
tributed amongst four p.^ orbitals, and the cluster orbital 
(with symmetric partners Vj-z^Vyz) is 




This cluster orbital aspect, adopted also by Erickson and 
collaborators, [8] is central to the following analysis. The 
overlap of these cluster orbitals is what is described by 

t(T : t^ , tg . 

Tendency toward magnetism. Xiang and Whangbo re- 
ported [18] obtaining a FM ground state for BNOO. 
Checking out this question with progressively finer fc-point 
samplings of the zone, we find BNOO is not quite unstable 
to a FM state. For meshes of 8x8x8 or finer, the mag- 
netic moment that is obtained on coarser meshes vanishes. 
Thus the Stoner interaction constant can be obtained from 
the condition of (near) fulfillment of the Stoner instabil- 
ity: 1st = 1/N{Ef) =0.35 eV [The Fermi level density 
of states N(Ei?)=2.88 states/eV per spin]. As will also 
be the case for the SOC strength (below), it is important 
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Fig. 3: (color online) Total and atom-projected DOS per for- 
mula unit, with spin-orbit coupling included. Oxygen contri- 
bution is 50% throughout the bands. The two lower bands 
are distinct, and nearly disjoint, leading roughly to a half- filled 
single-band configuration. 

that this parameter is not evaluated from the Os atomic 
potential alone, because the relevant local (molecular) or- 
bital is only half Os 5d, and 1st is reduced accordingly. 
The Hund's (exchange) splitting for a fully polarized 
orbital is Jh ~ 1st = 0.35 eV, which is smaller than the 
bandwidth but comparable to the SOC strength (see be- 
low). 

Spin-orbit coupling. It has been known for fifty years 
[19-21] that the t2g complex can be mapped onto an an- 
gular momentum L — > — £ with quantum number C=l, so 
within these bands the SOC operator becomes —£^£-S with 
^ > 0. The basis {Vxy^Vyz, V^x} can be mapped onto the 
Cz eigenstates to |0) — > V^y] | ± 1) ^ Vzx ± iT^yz, where 
the integer denotes the angular momentum projection mc 
onto z. The negative sign leads to an inverted spectrum, 
with the = I quartet at -^/2 and the J = \ doublet 
at -1-^. The calculated splitting at fc=0 gives a large SOC 
strength ^=0.30 eV, so large it splits off the entire upper 
J =1/2 band as is clear from fig. 1. 

The DOS with SOC included is displayed in fig. 3. In 
addition to the split-off upper J = \ band, the lower two 
bands are nearly disjoint, being bound together only by 
the degeneracy within the J = ^ states at zero momen- 
tum(fig. 1). As a result, the problem of the Mott insulat- 
ing state does not rigorously reduce to a single band prob- 
lem though it may be effectively single band. As a result 
of the mapping L — > —C within the t2g states of interest, 
the total angular momentum is, J = S — C {J = 5,5), the 
magnetic moment is M = /is (25* — £), and the factor is 
(7,/ = I — sj(J+i) ■ ^^^^ isolated cluster limit, the r.m.s. 
moment (Af^)^/^ is small (1 /is) for J = ^ but is larger 
{^/Ihb) in the -/ = | ground state. 

At the band structure level, there remains only the 



Hund's exchange splitting H^x = ^JhSz due to the Os 
moment, which becomes important because this interac- 
tion lifts the remaining degeneracies. The lowest state at 
A:=0 is 

\mcms) = I + 1 T) « \[Dzx + iVyz] T). (2) 

but the dispersion leads to overlapping bands and metallic- 
ity, so intra-cluster Coulomb repulsion effects are required 
to describe the observed insulating character, as also con- 
cluded by Xiang and Whangbo. 

Intra-atomic versus intra- molecular repulsion: the Mott 
insulating state. For a lattice of d} ions in t2g orbitals, 
an on-site Hubbard repulsion of roughly U ^ W is ex- 
pected to lead to a Mott insulating state. To model this 
insulating state we initially applied the correlated band 
LDA-I-U method, using J=l eV (whose value seemed im- 
material) and varying U. Since the initial interest was to 
look for an orbitally ordered solution, the full rotationally- 
invariant form of interaction implemented in Wien2k was 
used, and (following conventional practice) the U (and J) 
interaction is applied only to the Os 5d orbitals. A Mott 
insulating state was not obtained for values of U even sev- 
eral times larger than W. At U =4 eV « 4W the increasing 
exchange splitting leads only to a half metallic state, not 
an insulating state. 

This failure of the LDA-I-U method, as conventionally 
applied to atomic orbitals, to produce the expected Mott 
insulating state is a situation we have observed previously 
in Ad and 5d oxides. The problem can be cast not as 
a failure of the LDA-f U approach per se, but rather of 
its application to an extended cluster orbital by applying 
the orbital potentials solely to the Os 5d atomic function, 
rather than the entire molecular orbital T>. Since only 50% 
of the V orbital's charge is Os 5d, the LDA-I-U correction 
tends mainly to change the relative amounts of 5d and 
2p character, rather than to open a gap by shifting the 
energy of a cluster orbital. As a result, no state with 
integer 5d occupation number can arise, whereas in the 
LDA+U method integer occupation is necessary to obtain 
the insulating state. This problem is cured by using the 
Wannier basis, rather than the atomic state basis, in the 
LDA+U method. 

The full correlated electronic structure. The unfolding of 
the electronic structure that we find most instructive is 
pictured in fig. 4. We begin from the basis of degenerate 
orthonormal cluster orbitals (Hq), include SOC interac- 
tion Hsoc, then account for Hund's (exchange) interac- 
tion that breaks the degeneracy, and finally the repul- 
sive Coulomb interaction Hjj (as in LDA-I-U) 

H = Ho+^L-S~JhSz + Hu. (3) 

All of these are single site, and therefore simple to visual- 
ize. Kinetic effects (banding) is a small effect that enables 
charge transfer, but when it does not destroy the gap, as 
here, it is the least significant. 
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Table 1: Values of spin, orbital, and total moments versus the 
value of U , within the correlated t2g shell of BNOO. Increasing 

the hopping t parameters (by 50-100%) reduces the individual 
moments but hardly changes the net moment. 



u 


Ms 


Ml 


Mtot 


1.1 


0.714 


-0.717 


-0.003 


1.5 


0.909 


-0.866 


-1-0.043 


2.0 


0.976 


-0.935 


-1-0.041 



As mentioned above, Hund's exchange lowers the pure 
spin up I -I- 1 t) niember of the lower spin-orbit quartet, 
making it the orbital to occupy to achieve the lowest en- 
ergy. Then Coulomb repulsion lowers this occupied state 
by ~ [/ with respect to the others. We find that a critical 
(minimum) value olUc — 1-1 eV is required to produce the 
Mott insulating state. This critical value is significantly 
less than (therefore consistent with) values that have been 
suggested [8,18] as applicable to BNOO. This agreement 
reafiirms that Coulomb repulsion on top of strong SOC 
is an essential aspect of the mechanism underlying the 
insulating state. The compensation of the moments dis- 
tinguishes between the various viewpoints, as we now de- 
scribe. 

The I + 1 t) state is pure-spin and pure-orbital mo- 
ment, and within this state the z-component of the mo- 
ment (M^) = (-1-1 t \{2Sz-C^)\ + l t) = is compensated 
precisely by opposing spin and orbital contributions. This 
feature makes BNOO an unprecedented transition metal 
analog of the case of the Sm'^^ c?^ ion [22] where L = 2S 
and Mz ~ 0: the ion/cluster is magnetic (J 7^ 0) but 
there is nearly vanishing total moment for a fixed direc- 
tion of spin. The effect of mixing-in of other states by the 
kinetic hopping (quenching) is measured by the resulting 
exchange coupling 4:jj ~ 30 meV {t ~ 120 meV, U=1.5- 
2 eV), i.e. it is small compared to other energies in the 
system. While the spin and orbital moments each depend 
fairly strongly on the value of U (table 1), the cancella- 
tion between the two remains near-complete. This cancel- 
lation is a nontrivial effect, since appreciable mixing-in of 
the )0 I) state along (say) would degrade the spin moment 
without affecting the orbital moment. 

Discussion. The small ordered moment was one of the 

principle questions about BNOO. From the point of view 
presented here, a small moment is expected and the ques- 
tion becomes rather: why is it as large as observed? 
Naively, this near-pure spin state is unexpected, given the 
strong SOC which is largely regarded as mixing {not sepa- 
rating out) spin and orbital moments, nevertheless it is an 
occurrence that has recently been anticipated. [23] At the 
level of description of the Hamiltonian eq. 3 in mean field 
approximation, we have a Curie- Weiss moment (M^)^/^ 
= \/7M-B and an ordered FM moment of ss 0.04/iB (see 
table 4), which is only nonzero due to quenching. The 
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Fig. 4: Evolution (schematic) of the t2g energies from isolated 
3 X 2-fold degenerate for the isolated ion (in its cubic environ- 
ment), through SOC, Hund's coupling. Coulomb repulsion U , 
and finally dispersion due to hopping t. The occupied band has 
pure spin-orbital | -|- 1 |) character until intercluster hopping 
is included. 

agreement with the observed values (near 0.7 and 0.2 ^b, 

respectively) in absolute terms is not good, however one 
would expect a small ordered moment of 0.2 /xb to be diffi- 
cult to reproduce in a first principles manner. The degree 
of coUinearity of S and C, which finally determines the 
moment, is sensitive to the environment. The magnitude 
of collinearity, || £ • S || /(|| 5 |||| £ ||), is only 40% in the 
J = I manifold anyway. Adding a crystalline anisotropy 
term decreases it rapidly, by a further 35% for 7 = ^ 
and by 60% for 7 = 2^ (this is orbital moment quenching). 
The quantum fluctuations are thus quite sensitive to the 
environment. 

The matter of quantitative agreement of moments 

is actually a broader question: the related compound 
Ba2CaOs06 (BCOO) with one more electron and 5=1 and 
a cubic crystal structure is more conventional, yet it has a 
Curie- Weiss moment [24] of 1.61 /is that is also far below 
the spin-only value of 2.83 jiB- The Re-based d} dou- 
ble perovskites Sr2XRe06, X = Mg and Ca, on the other 
hand, experience much smaller SOC and have moments 
equal to the spin only values and display structural dis- 
tortions as the conventional picture would suggest. This 
distinction supports the involvement of strong spin-orbit 
coupling in both the small moment and the "restoration" 
of cubic symmetry in the heptavalent osmates BNOO and 
BLOO. 

The magnetic ordering issue (BNOO is FM; isostruc- 
tural and isovalent BLOO with almost identical volume 
is AFM) is a challenging and delicate question involving 
several complications. First, the fee lattice itself can be 
considered one of the most frustrated. Frustration on the 
triangular lattice, with its edge-sharing triangles, is well 
known. Symmetric tetrahedra are frustrated, and the py- 
rochlore lattice of vertex-sharing tetrahedra is a well stud- 
ied case of strong frustration. The fee lattice can be pic- 
tured as one of space-filling, face-sharing tetrahedra. 
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Second, spin coupling between the OsOe clusters in 
BNOO proceeds through three distinct hopping processes 
ta,t„,ts between anisotropic V^z + i'Dyz cluster orbitals, 
and the behavior of the Curie- Weiss susceptibility re- 
flects a net antifcrromagnctic coupling. Such coupling is 
strongly frustrated on an fee lattice, and the actual pattern 
of ordering can be a delicate problem. The difference be- 
tween FAf BNOO and AFM BLOO, based on some small 
detail, reflects the geometrical frustration and apparent 
competing couplings that must be expected in this struc- 
ture. Note that the cousin S=^ compound Sr2CaRe06, 
even with the coherent structural distortion that relieves 
some of the frustration, has been shown to be a spin glass 
system. [26] 

Thirdly, orbital ordering is also subject to frustra- 
tion. [27] The surviving cubic structural symmetry of the 
Os ion is seems remarkable, but one must recognize 
that the "ferromagnetism" arising in a strongly coupled 
spin-|-orbital moment system lowers the electronic sym- 
metry of the system. The question of whether this degree 
of freedom remains disordered (a combined 'spin-orbital 
liquid' within the FM phase), and thus serves to restore 
effective cubic symmetry, is a question that requires fur- 
ther study. In any case, frustration and large spin-orbital 
coupling are key issues in Ba2NaOs06. 

This work was supported by DOE grant No. DE-FG03- 
01ER45876 and DOE's Computational Materials Science 
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